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ABSTRACT 



Context. Mid-infrared (mid-IR) emission lines of molecular hydrogen (Ht) are useful probes to determine the mass of warm gas 
present in the surface layers of circumstellar disks. In the past years, numerous observations of Herbig Ae/Be stars (HAeBes) have 
been performed, but only two detections of H2 mid-IR emission toward HD 97048 and AB Aur have been reported. 
Aims. We aim at tracing the warm gas in the circumstellar environment of five additional HAeBes with gas-rich environments and/or 
physical characteristics close to those of AB Aur and/or HD 97048, to discuss whether the detections toward these two objects are 
. suggestive of peculiar conditions for the observed gas. 

' Methods. We search for the H2 S(l) emission line at 17.035 /jm using high-resolution mid-IR spectra obtained with YLI/VISIR, 

and complemented by CH molecule observations with VLJ/UVES. We gather the H2 measurements from the literature to put the new 
, results in context and search for a correlation with some disk properties. 

^ |i Results. None of the five VISIR targets shows evidence for H2 emission at 17.035 fim. From the 3cr upper limits on the integrated 

I ' line fluxes we constrain the amount of optically thin warm (> 150 K) gas to be less than ~ 1.4 Mjup in the disk surface layers. There 

O ' are now 20 HAeBes observed with VISIR and TEXES instruments to search for warm H2, but only two detections (HD 97048 and 

AB Aur) were made so far. We find that the two stars with detected warm H2 show at the same time high 30/13 yum flux ratios and 
C/j . large PAH line fluxes at 8.6 and 1 1.3 fim compared to the bulk of observed HAeBes and have emission CO lines detected at 4.7 fjm. 

^ ' We detect the CH 4300.3 A absorption line toward both HD 97048 and AB Aur with UVES. The CH to H2 abundance ratios that this 

would imply if it were to arise from the same component as well as the radial velocity of the CH lines both suggest that CH arises 
from a surrounding envelope, while the detected H2 would reside in the disk. 

Conclusions. The two detections of the S(l) line in the disks of HD 97048 and AB Aur suggest either peculiar physical conditions 
or a particular stage of evolution. New instruments such as Herschel / PACS should bring significant new data for the constraints of 
\^ \ thermodynamics in young disks by observing the gas and the dust simultaneously. 

I Key words, stars: circumstellar matter - stars: formation - stars: pre-main sequence ~ ISM: molecules 

1. Introduction Therefore, molecular line emission is not sensitive to and does 

not probe the mid-plane interior layers of the disk because it is 

Planets are supposed to form m circumstellar disks composed of ^ j^^y ^^^^y^ g^^^^^^ ^^^^^^ molecular gas in the opti- 

gas and dust around stars in their pre-main sequence phase. At ^^jj^ ^^^^ ^^^^^^^ j^y^^.^ ^^^^jj^ t he expected H, line fluxes ai-e 

this evolutionary stage, the disk mass is essentially dominated by ^^^j^ ^„ ^^^^^pj^ ^^^^^ ICarrnona et alJ (l2008h calcu- 

the gas (99%), especially molecular hydrogen (H2). Although H2 j^jg^ j^e expected H2 line fluxes from typical Herbig Ae disks 

IS the principal gaseous constituent in disks, it is very challeng- ^^^^^^^^ ^ two-layer Chiane & Goldreich ( 1997) disk model, 

ing to detect. Molecular hydrogen is an homonuclear molecule j^^^^ ^ j^^^^^ ^ gas-to-dust-ratio of 100, LTE emission and a dis- 

which means that its fundamental transitions are quadrupolar in ^^^^^ ^1^^^ expected line H2 fluxes are 

nature. Hence their Einstein spontaneous emission coefficients ^^^^ ^^^^^^^ ^^^^ detection limits of cuiTent instrumentation, 

are small and produce only weak lines. For circumstellar disks j^^^g^^ numerous non-detections of H2 mid-IR pure rotational 

another chaUenge is that the weak H2 hnes must be detected on jj^^^ j„ circumstellar environment of young stars have been 

top of the strong dust continuum emission. In the mid-infrai-ed, reported in flie past few y ears (Bitner et al. 2008; Ca rmona et all 

high specti-al resolution instruments like VISIR at the VLT ai-e ^ iMartin-ZaMi et al.ll2008k ,2009 a). Nevertheless, lines can 

required to disentangle these weak lines from the continuum. ^^^^j^ ^ detectable level if flie gas-to-dust to ratio is allowed to 

However, one should bear in mind a fundamental issue con- ^e higher than 100 and/or the T„,, > Tdus, in the surface layers 

cerning the structure of a gas-rich optically thick disk. Molecular j-j^g ^j^j^ 
lines are produced in the hot upper surfaces of the disk where 

the molecular gas and accompanying dust are opticaUy thin. Molecular hydrogen mid-IR lines have been detected in two 

Herbig Ae/Be stars, namely HD 97048 and AB Aur, respectively 

Send offprint requests to: C. Martin-Zaidi observed with VLT/VISIR (.Martin-Zaidi et alj2007h and TEXES 
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Table 1. Astrophysical parameters of the sample stars. 



Star 


Sp. 


<■// 




(a) 


J 

a 


Disk 


nviuence tor 




Type 


(K) 


(mag) 


(km s-') 


(pc) 


resolved 


gas-rich 
















environment 


HD 142527 


F6 Ille 


6300") 


1.49"' 


-3.5*-^' 


140(" 


yes'"*' 




HD 169142 


A8 Ve 


8130'^' 


0.37<« 


-3.0<'' 


145<5> 


yes'*' 


yes***' 


HD 150193A 


Al Ve 


9300<" 


1.61"' 


-6.0"'» 


150"' 


yes"" 




HD 163296 


Al Ve 


9300"' 


0.25"' 


-4.0<2' 


122"' 


yes"2' 


yes"3' 


HD 100546 


B9 Vne 


10470"' 


0.25"' 


+ 17(13) 


103"' 


yes(i4) 


yes"3> 



'"' radial veloc ity of the star in the heliocen tric rest frame, 

References: "' Ivan den Ancker et al.llll998h: <^' S IMBAD d atabase: '^' [ de Zeeuw et an(ll999h 

j200l[) : '^' iPanic et al.' (f200l 



<^' iMalfait et all (Il998h : iDunkin et all (|1997|): ''"iKuhn et al 
(I2003h : "^' iGradv et alJ (l2000h : "^' iLecavelier des Etangs et al 



20031) : "*'|Augereau et al._ 



<^'|Fukagawa et al. (2006^; '^' lAcke et al.](l2005h: 
^; "" ' .Reipurth et ah U99l) : "" iFukagawa et all 
200 Ih. 



Table 2. Summary of the observations. 



Star 


^exp 


Airmass 


Optical 


Standard 


Airmass 


Optical 


Asteroid 


Airmass 


Optical 




(s) 




Seeing 


Star 




Seeing 






Seeing 








(") 






(") 






(") 


HD 142527 


3600 


1.10-1.28 


0.71-1.08 


HD 21 1416 


1.23-1.25 


0.79-1.02 


HEBE 


1.19-1.20 


0.97-1.06 


HD 169142 


3600 


1.17-1.50 


0.71-1.01 


HD 21 1416 


1.23-1.25 


0.79-1.02 


HEBE 


1.19-1.20 


0.97-1.06 


HD 150193A 


3600 


1.31-1.88 


0.85-1.40 


HD 21 1416 


1.24-1.25 


0.81-0.99 


PSYCHE 


1.13-1.16 


0.75-0.83 


HD 163296 


3600 


1.60-2.70 


0.79-1.19 


HD 211416 


1.24-1.25 


0.81-0.99 


IRIS 


1.00-1.01 


0.99-1.28 


HD 100546 


3600 


1.44-1.52 


0.73-1.09 


HD 89388 


1.31-1.32 


0.61-0.73 


HEBE 


1.19-1.20 


0.97-1.06 



The airmass and seeing intervals are given from the beginning to the end of the observations. 



dBitner et alJl2007l) . These detections imply particular physical 
conditions for the gas and dust, such as Tg^s > Tdust, as men- 
tioned above. These conditions may be created by gas heated 
by X-rays or UV photons in the surface layers of the disks. We 
note that these two detections indicate that the gas has not com- 
pletely dissipated in the inner part of these disks in a lifetime of 
about 3 Myrs (ages of the stars), while photoevaporation of the 
gas is expected to clear up this inner region wi thin roughly th e 
same ti me (<3 Myrs; e. g. Takeuchi et al. 2005; Alexanderl2008h . 
Indeed, Martin-Zaidi e t al.i ((2009a) have shown that the emitting 
H2 around HD 97048 was more likely distributed in an extended 
region wi thin the inner disk , between 5 AU and 35 AU of the 
disk, and iBitner et al ] (120071) concluded from their observation 
of the disk of AB Aur that the H2 emission arised around 18 AU 
from the central star. 

We present new observations with the VISIR high-resolution 
spectroscopic mode, to search for the mid-IR H2 emission line 
at 17.035 fim, the most intense pure rotational line observable 
from the ground, in a sample of well studied Herbig Ae/Be stars 
known to harbor extended gas-rich circumstellar disks. The main 
goal of these observations is to enlarge the global sample of 
HAeBes observed at 17.035 fim, and to better constrain the par- 
ticular physical conditions observed in the disks of HD 97048 
and AB Aur In order to better understand the detections of H2 in 
the disks of HD 97048 and AB Aur, we also analyze VLT/UVES 
spectra of these two stars to observe spectral lines of the CH 
and CH^ molecules which are linked to the formation and ex- 
citation of H2. In Sects. |3]and|4]we present the analysis of the 
VISIR spectra of the five HAeBes of our sample and relate it to 
the context of the global search for molecular hydrogen at mid- 
IR wavelengths. We perfom a statistical analysis of the whole 
sample of HAeBes where H2 emission has been searched for in 



the mid-IR, and explore the possible link between the H2 and the 
dusty disk properties. Finally, we discuss the possible origin of 
the mid-IR H2 emission in the disks of HD 97048 and AB Aur. 

2. Observation and data reduction 

2.7. VISIR observations 

We selected a sample of five Herbig Ae/Be stars with spatially 
resolved disks (Table [T]). They were selected using the follow- 
ing criteria: (i) proximity to the Sun (<150pc), with well known 
distances; (ii) resolved disk and/or firm evidence of gas-rich 
circumstellar environment, and/or properties close to that of 
HD 97048 and/or AB Aur, for which the 17.035 //m line has al- 
ready been detected; ( Hi) well known Herbig Ae/Be stars with a 
rich set of existing observations, particularly including gas-lines 
spectroscopy ob s ervations (Far-UV H2 line s , CO lines...; e.g. 
Thi et all l200U iMartin-Zaidi et alJ l2008at Ivan der Plas et al.1 
20091) and R\Hs (Polycy c hc Aromatic Hydrocarbons) fe atures 
(e.g. iHabart et alj 12004 lAcke & van den Anckeii l2006l) . The 
fundamental parameters of each target are presented in Table [T] 
The target stars were observed on 2009 June 14 and 15, 
with the high spectral resolution long-slit mode of the VISIR 
(ESQ VLT Image r and Spectrometer for the mid-InfraRed 
\Laeaee et al. 12004) instrument at the VLT. We focussed on the 
observations of the S(l) pure rotational line of H2 at 17.035 //m 
because it is the most intense line observable from the ground, 
and because it is the onl y line ever detected with V ISIR in disks 
of Herbig Ae/Be stars jMartin-Zaidi et al.1 12007|) . The central 
wavelength of the observations was thus set to be 17.035 jum. 
We used the 0.75" slit, providing a spectral resolution of about 
14000. The observation conditions are summarized in Table |2] 
In order to correct the spectra from the Earth's atmospheric ab- 
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- obs- HD142527 

- Obs. HEBE 



1 7,02 1 7.03 1 7.04 
Wovelenglh {fim) 




1 7.02 1 7.0J 1 7.04 
Wovelenglh ^tm) 




1 7.0J2 1 7.034 1 7.036 1 7.038 1 7.040 
Woveiength (jim) 




- obs. Hni691-1? 

- Obs. HEBE 




iZ 6 
5 





1 3a 











1 7.02 1 7.03 1 7.04 
Wovelenglh (/Am) 



1 7.02 1 7,03 1 7.04 
Wovelenglh i/tm) 



17.032 17.034 17.036 
Woveiength {/tm) 






i 
i 
; 
: 
: 




■ 1 



n ; 3a 










I 


■ ■ 1 



1 7.02 1 7.03 1 7.04 
Wovelenglh (/im) 



1 7,02 1 7.03 1 7.04 
Wovelenglh (ftm) 



1 7.032 1 7.034 1 7.036 1 7.038 1 7.040 
Woveiength (ftm) 




obs. HDl 63296 
Obs. IRIS 





1 7.02 1 7,03 1 7,04 
Wovelenglh (firn) 



1 7,02 1 7,03 1 7.04 
Wovelenglh (/tm) 



17,032 17,034 17.036 17,038 17.040 
Woveiength {nm) 




1 1 1 


1' ' ' 1 1 ' ' ', 








1 7.02 1 7.03 1 7.04 
Wovelenglh (^m) 



1 7.02 1 7.03 1 7.04 
Wovelenglh {ftm) 



17.032 17,034 17.036 17.038 17.040 
Woveiength (ftm) 



Fig. 1. Spectra obtained for the H2 S(l) line at 17.035 //m. Left panel: continuum spectra of the asteroid and of the target before 
telluric correction. Middle panel: full corrected spectra: dotted lines show spectral regions strongly affected by telluric features. 
Right panel: zoom of the region where the H2 Unes should be observed (dashed vertical Unes). A Gaussian of width FWHM = 21 
km s"^ and integrated line flux equal to the Scr line-flux upper limits is overplotted. The spectra were corrected neither for the radial 
velocity of the targets nor the Earth's rotation velocity. 
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sorption and to obtain the absolute flux calibration, asteroids and 
standard stars were observed at nearly the same airmass and 
seeing conditions as the objects, just before and after observ- 
ing the science targets. For this p urpose, we also use d the mod- 
eled spectra of the standard stars (ICohen et alJll99S j), and those 
of the asteroids (for HEBE and IRIS see Mueller & LagerrosI 
Il998ll2002l and for PSYCHE, Muefler 2009, priv. comm.). For 
details on the observatio ns a nd data reduction techniqu es see 
iMartin-Zaidi et al.1 ( l2007h and lMartin-Zaidi et all (l2008bl) . 



2.2. UVES obseryaimns 

To better understand the detections of H2 towards HD 97048 
and AB Aur, we observed spectral lines of the CH and CH^ 
molecules in the optical range at high spectral resolution using 
WES (ESO program numbers 075.C-0637 and 078.C-0774). 
The blue arm of UYES (3700A - 5000A) gives acces to the 
electronic transitions of the CH and CH^ molecules that are 
good tracers of the H2 formation an d excitation ( Federman 1982 
lMattila"l986- Somervifle & Smith T989). The formation of CH 
is predicted to be controlled by gas-phase reactions with H2. 
CH is thus a good tracer of H2 and their abundances are gen- 
erally strongly correlated. The formation of the CH^ molecule 
through the chemical reaction C^ + H2 needs a temperature of 
about 4500 K to occur. Thus, the CH^ molecule is a probe of 
hot and excited media, which could be interpreted as material 
close to the star, and this would allow us to better constrain the 
excitation of H2. 

Th e sp ectra were reduced using the UVES pipeline v3.2.0 
(iBallester et al.l I2OOOI) based on the ESO Common Library 
Pipeline, and were corrected from the Earth's rotation in order to 
shift them in the heliocentric rest frame. To analyze the absorp- 
tion lines in the observed spectra we used the Owens package 
(for details, see Lemoine et al. 2002; Martin-Zai'di et al. 2008al), 
which allows us to derive the characteristics of the gaseous com- 
ponents, i.e. column density N, radial velocity Vrad and broaden- 
ing parameter b. 



3. Results 

3.1. VISIR sample analysis 

We observed five well studied Herbig Ae/Be stars with the high 
resolution spectroscopic mode of VISIR to search for the H2 pure 
rotational S(l) line. None of the observed sources show evidence 
for H2 emission at 17.035 yum (Fig.[B. In HD 169142, the ex- 
pected centroid for the H2 line corresponds to an emission fea- 
ture in the spectrum. However, this feature has a full-width-at- 
half maximum (FWHM) smaller than a spectral resolution ele- 
ment (Av ~ 2 1 km s"' ), and an amplitude of about 1 cr that makes 
it insignificant in terms of detection (for detail s on the calcula - 
tion of the standard deviation cr, see Martin-Za idi et al ] l2008bh . 
We did not consider it further as a detection. As described pre- 
viously, we insist that mid-IR H2 lines only probe warm gas 
located in the surface layers of the disks because of the large 
opacity in the interior layers due to the dust. The surface lay- 
ers of the disks, from which solid-state emission features also 
arise, surrounding our target stars therefore do not contain sufli- 
cient warm gas to enable it to be detected in emission at mid-IR 
wavelengths. 

Then, we derived 3cr upper limits on the integrated line 
fluxes and upper limits for the total column densities and masses 
of H2 as a function of the temperature for all the spectra by in- 
tegrating over a Gaussian of FWHM corresponding to the VISIR 



spectral resolution, and an amplitude of 3cr, centered on the ex- 
pected wavelength for the S(l) line (see Tableland Fig.|2l). 

For this purp ose, we used the method detailed in 
Ma rtin-Zaidi et al.l (2008b, 2009a), by assuming that the line is 
optically thin at local thermodynamic equilibrium (LTE) and that 
the radiation is isotropic (for d etails on the formulae, see also 
Ivan Dishoecklll992HTakahashill2001l) . 



100.000 



10.000 



1.000 




HD142527 

HD169142 

HD150193A 

H0 163296 

HO 100546 



600 800 

V (K) 



Fig. 2. Upper mass limits of optically thin H2 derived from H2 
S(l) line flux limits as a function of the assumed LTE tempera- 
ture. 



From the 3cr upper limits to the emission line flux, we cal- 
culated upper limits on the column density and mass of H2 for 
each star. We found that the column densities should be lower 
than -lO-"* cm^ at 150 K, and lower than -lO^^ cm^ at 1000 K. 
The corresponding upper limits to the masses of warm gas in 
the surface layers of the inner disk were estimated to be in the 
range from ~6xl0~^ to ~1.4 Mjup (1 Mjup~10"-' Mq), assuming 
LTE excitation, and depending on the adopted temperature (see 
Tableland Fig.©. 



3.2. H2 detection statistics 

In the past three years, numerous observations have been 
performed to observe the pure rotational mid-IR emission 
lines of molecular hydrogen in the circumstellar environ- 
ments of HAeB es. VISIR o bservatio ns of 15 HAeBes have 
been conducted ( Carmona et al. 200S: iMartin-Zaidi etal ] l2007t 
iMartin-Zaidi et alT 2008bl and this work) and only one source 
out of the 15, namely HD 97048, presents a clear evidence for 
H2 emission at 17.035 fim. In addition, in t heir sample of five 
HAeBes observed with TEXES at 17.035 ^m. lBitner etan(l2008h 
have reported only one detection of the S(l) H2 line in the disk 
of AB Aur. This leads to a detection statistics of the S(l) line 
from the ground of 10% only in disks about Herbig Ae/Be stars. 
However, this statistic suffers from the fact that the 3cr limits 
in our sample are inhomogeneous. Indeed, due to the differ- 
ent detection limits of each instrument (and exposure time for 
each target), the 3cr upper limits on the line flux from VISIR ob- 
servations cannot be directly compared to those obtained with 
TEXES. However, our sample is quite representative of the vari- 
ety of HAeBes observable from the ground, including stars with 
spectral type from F6 to B2, with very young, massives or tran- 
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Table 3. Results of the H2 S(l) line analysis. 



Star 
HD 


obs 

(jum) 


Integrated flux"'' 
(ergs s"' cm"^) 


Intensity"'' 
(ergs s"' cm"- sr"') 


A'(H2) upper limits*' 
{cmr-) 


TT 

"2 


mass upper limits"' '^' 
(Mj„p) 


150K 


300K 


lOOOK 


150 K 


300 K 


1000 K 


142527 


17.0351 


<1.0xl0-'* 


<2.5xlO-^ 


9.4x10-^ 


6.1x10^' 


1.4x10-' 


1.9x10- 


' 1.2x10-^ 


3.3x10-3 


169142 


17.0343 


<1.8xl0-''' 


<4.2xl0-3 


1.6x10-^ 


1.0x10^- 


2.4x10-' 


3.5x10- 


' 2.1x10-^ 


6.1x10-3 


150193A 


17.0349 


<2.1xl0-"' 


<5.0xlO--'' 


1.9x10" 


1.2x10^- 


2.8x10-' 


4.5x10- 


' 2.7x10-2 


7.7x10-3 


163296 


17.0345 


<2.8xlO-i'' 


<6.7xl0-3 


2.5x10^3 


1.6x10^2 


3.8x10-' 


3.9x10- 


' 2.4x10-2 


6.9x10-3 


100546 


17.0363 


<1.4xl0-'' 


<3.3xlO-^ 


1.2x10^^ 


7.9x10^- 


1.9x10-^ 


1.4 


8.4x10-2 


2.4x10-2 



(a) expected position of the line in the observed spectra by correcting the radial velocity of the star for the epoch of observation. 

(b) 3cr upper confidence limits. 

(c) masses of H2 are calculated assuming the distances quoted in Table[T] 



sitional disks. Therefore, the 10% statistics can be regarded as a 
lower limit on the number of H2 detections with current ground- 
based mid-IR instrum e ntatio n. This is consistent with the mod- 
els by ICarmona et al. I (l2008l) . who show that mid-IR H2 lines 
cannot be detected with the existing instruments for disks under 
LTE conditions. Those authors estimated that sensitivities down 
to IQ-'^ ergss-' cm-^ must be reached to detect the S(l) line, 
while VISIR typically reaches IQ-'"* ergs s-' cm-^. 

The numerous non-detections of the S(l) line in circumstel- 
lar environments of HAeBes could imply that the physical con- 
ditions of the warm gas in most disks are consistent with those 
assumed in such a model, i.e. gas and dust well-mixed, a gas-to- 
dust ratio of about 100, and T„as - Tdusi- 



4. Discussion 

4.1. H2 and disk properties 

The challenge now is to explain the two detections of the S(l) 
line in the disks of HD 97048 and AB Aur. If the HD 97048 
and AB Aur disks trace a particular evolutionary stage, one may 
expect some connection between the mass of warm H2 and the 
disk properties. 

The slope of the spectra corresponding to the ratio of the F30 
and F12 indexes, i.e. c ontinuum fluxes at 3 and 13 /im respec- 
tively, has been used bv'Bro wn et al.l (|2007|) to identify the popu- 
lation of cold disks. In addition, young flared disks tend to show 
rising spectra in the mid-IR with a much higher F30 / F13 ratio 
than the more evolved, settled sys tems with self-shado wed disks 
with much flatter mid-IR spectra (iDominik et al.ll2003h . The F30 
/ Fi3 ratio diagnostic is therefore not unique, but allows us to 
identify a population with properties that depart from those of 
most young stars with d isks. Indeed, in their sample of T Tauri 
stars, iBrown et al.l (1200 7) have shown that the majority of their 
sample had an emission that increased by a factor of 2.3+1 .4 be- 
tween 13 and 30 //m, while the cold, possibly transitional, disks 
rose by factors of 5 - 15. 

Here we measured the slope of the spectra of our targets, ob- 
tained with the IRS spectrograph installed onboard the Spitzer 
Space Telescope, by computing F13 which is the mean flux value 
(in Jy) in the range 13 yum + 0.5 /vm, and F^o the mean flux 
val ue (in Jy) between 29 fim and 3 1 fim (Table |4] for details 
see lOlofsson et all l2009l) . The data r eduction was perfo rmed 
with the "c2d legacy team pipeline" (iLahuis et al ] l2006h with 
the S18.7.0 pre-reduced (BCD) data. 

Figure [3] shows the upper limits of H2 mass measured from 
the S(l) line flux versus the slope of the spectra (Fso/Fb ratio). 



Table 4. ^30 / F13 indexes, and H2 mass upper limits at 600 K. 



Star 


^30 / -^^13 


H2 mass and upper 






limits at T=600 K (Mj„p) 


(This work) 






HD 142527 


4.86 


<4. 1x10-3 


HD 169142 


7.12 


<7.5xlO-3 


HD 150193A 


1.44 


<9.6xl0-3 


HD 163296 


1.54 


<8.5xlO-3 


HD 100546 


3.32 


<2.9xl0-2 


(Martin-Zaidi et al. 2008b) 






HD 98922 


0.70 


<7. 1x10-2 


HD 250550 


2.67 


<1.2xl0-' 


HD 259431 


2.78 


<3. 1x10-2 


HD 45677 


(a) 




(Carmona et al. 2008) 






UX Ori 


(a) 




HD 100453 


5.41 


<3. 8x10-3 


HD 101412 


0.92 


<7. 3x10-3 


HD 104237 


1.28 


<3.5xlO-3 


HD 142666 


1.55 


<4.7xl0-3 


(Bitner et al. 2008) 






49Cet 


2.31 


<6.6xl0-'' 


MWC 758 


4.00 


<5. 3x10-3 


V892 Tau 


3.72 


< 1.9x10-2 


VVSer 


0.76 


<1.5xl0-2 


(Martin-Zaidi et al. 2009a) 






HD 97048 *' 


6.22 


(1.6±0.1)xl0-2 


(Bitner et al. 2007) 






AB Aur "'' 


4.21 


(1.6±0.47)xl0-3 



(a) Spitzer data not available. 

(b) masses of H2 (not upper limits) because the S(l) has been de- 
tected. 

We used distances and integrated fluxes values from the papers ref- 
erenced in the table. 



For comparison with our targets, we plotted all the stars observed 
at 17.035 pm by VISIR and by TEXES, including HD 97048 and 
AB Aur (Martin-Zaidi etal. 2007; Martin-Zaidi et al. 2008^ 
l2009al: ICarmona et al.1 120081: iBitner et al.l i2008). For this task, 
we considered only the upper limits on the mass of H2 at 600 K 
for each target, because the temperature of the observed H2 is 
around 600 K for AB Aur (Bitner et al. 2007) and lower than 
600 K for HD 97048 (Martin-Zaidi et al. 2009a). 

Most of our sample stars have an emission that increases by 
a median factor of 2.2+1.5 between 13 and 30 pm, which is 
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Fig. 3. Slope of the Spitzer spectra (F20/F12 ratio) versus upper 
mass limits of optically thin H2 at 600 K derived from H2 S(l) 
line flux limits. In this diagram are plotted all the Herbig Ae/Be 
stars observed at 17.035 fim with TEXES and VISIR. Triangles 
represent the two stars for which the S(l) line has been detected. 
Squares represent the Herbig Be stars for which a circumstellar 
disk has never been clearly detected. Plain circles represent the 
Ae stars of the sample (see text). 



fully consistent with the results of lBrown et al ] (l2007h . However, 
four stars in our sample, namely HD 169142, HD 100453, 
HD 142527, and HD 97048, have significantly higher Fjo/Fn 
ratios. We stress that AB Aur is marginally higher than the bulk. 
In addition, the four spectra are also characterized by the ab- 
sence of 10 fj.m amorphous silicate features. This is consistent 
with the fact that the strong wings of the amorphous silicate 
feature at 10 yum increase the continuum flux around 13/vm and 
thus decrease the F^o/Fn ratio. The F^o/Fn ratio of these four 
stars c orrespond to the region of the plot where iBrown et alJ 
(I2OO7I) have identified the transitional cold disks in their sample, 
which also do not present the 10 jum amorphous silicate feature 
in their spectra. However, HD 97048 is sup posed to be a young 
star harboring a ma ssive gas-rich flared disk ('Lagag e et al.l2006t 
iDoucet et al.l 12007^ . as compared to the other three stars that 
are known to have more evolved and less massive disks than 
HD 97048 (e.g. Grady et al. 2007). Its position in our Fig. [3] 
raises numerous questions about its status. Below, we explore 
different ways to understand the origin of the detected H2 and 
the status of HD 97048 (and marginally AB Aur). 



4.2. Is the H2 related to PAHs ? 

As we did for the Fso/Fn ratio, we searched for the possible 
correlation between the presence and line flux of PAHs features 
and the upper limits on H2 mid-IR lines. For this task we com- 
pared the line flux of the PAHs feature at 8.6jum and I1.3yum 
(called 1 lyum complex due to the possi ble blend with features of 
amorp hous and crystalline silicate; e.g. lAcke & van den Anckerl 
I2OO4I) found in the literature (Table |5]l, with the H2 upper mass 
limits (Table m. With regard to the F30/F13 ratio, AB Aur and 
HD 97048 seem to depart from the bulk, and appear in the top of 
the plots where the PAHs line fluxes are high (Fig.|4]and Fig.|5]). 

In the upper disk surface layers, photoelectric heating is very 
eflicient on small grains such as PAHs, and can play a signifi- 



Fig. 4. Line flux of the PAHs feature at 8.6yum versus upper mass 
limits of optically thin H2 at 600 K derived from H2 S(l) line flux 
limits. Same legend as for Fig.[3j 
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Fig. 5. Line flux of the PAHs complex at llyum versus upper 
mass limits of optically thin H2 at 600 K derived from H2 S(l) 
line flux Umits. Same legend as for Fig. |3] 



cant role in the gas heating process (" Kamp & Dunemondll2004t 
Jonkheid et al. 2007). The PAHs emission depends on the geom- 
etry of the disk, i.e. the emission is proportional to the amount 
of PAHs that is directly illuminated by the UV light from the 
star. Thus the PAHs emission depends both on PAHs abundance 
and on geometry, and it is unclear which of these two effects 
is responsible for the higher PAH emission strength in AB Aur 
and HD 97048 compared to other Herbig stars in the sample. 
Irrespective of its origin, the strong emission of PAHs appears 
related to H2 emission. Is there a link between their excitation 
and heating mecanisms ? This question must be studied further. 

As shown in Fig. |4] and Fig.|5] the spectra of the other stars 
of the sample present weaker PAH features. One could assume 
that for these stars photoelectric heating is less efficient and does 
not excite enough the H2 to be detected in the mid-IR range. 

In both plots HD 100546 also departs from the bulk. This 
may be because the disk of this star has a particular structure 
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Table 5. Line fluxes of the PAHs features at 8.6jum and line 



fluxes of the PAHs complex at 1 l/jm. 


star 


PAHs 8.6 (W/mO 


Comp 1 1 (W/m^) 


ref 


HD 142527 


<9.90xlO"'^ 


2.09x10-'" 


1 


HD 169142 


2.55x10"''* 


8.18x10-'^ 


1 


HD 150 193 A 


<3.64xl0"'' 


8.29x10-" 


1 


HD 163296 


<5.57xlO"'^ 


1.49x10-'^* 


1 


HD 100546 


5.61x10""* 


3.89x10-'^ 


1 


HD 98922 


y 


1.4x10-'"* 


2, 3 


HD 250550 


n 


n 


4 


HD 259431 


— 


- 


- 


HD 45677 


— 


— 


— 


UX Oi l 


<1.57xlO-'' 


<7.20xl0-'* 


1 


HU 1U045J 


<1. 65X11) 




1 
1 


HD 101412 


y 


1.5x10-" 


2,3 


HD 104237 


<2.67xl0"'^ 


<1.56xl0-" 


1 


HD 142666 


<5.37xlO-'^ 


5.42x10-" 


1 


49 Get 








MWC 758 








V892 Tau 








VV Ser 


<1.67xl0-'^ 


2.04x10-" 


1 


HD 97048 


4.52x10"''' 


3.97x10-''* 


1 


AB Aur 


3.68x10-"* 


4.61x10-"* 


1 



- Data not available. 

"Y" (yes): presence of PAHs features but flux not measured, 

"N" (no): absen ce of PAHs features. 

Refere nce s: 1: |Acke & van den Ancked d2004h: 2: iGeers et all 
( |20M) : 3: IGeers et alj j2007b : 4: lHabart et alj ( l20()4h ; 



with an inner ring, a gap between ~0.5 aiid 10 A U, and an outer 
disk (iBouwman et al.ll2003l:lBenistv et al.ll201(]h . 

4.3. Could the H2 emission arise from the envelope? 

Both HD 97048 and AB Aur are known to possess extended 
envelopes with a significant contribution in the various obser- 
vation s of gas and dust (e.g. iHartmann et al.llT993l: iGradv et"al] 
Could the observed H2 mid-IR emission be due to the 
gaseous component of the envelopes of these stars? 

This scenario can certainly excluded because, firstly, the H2 
S(l) fines observed in the spectra of HD 97048 and AB Aur 
are not spatially resolved, which constrains the emitting region 
to be very close to the central stars. Given the weakness of 
the detected H2 emission, at least for AB Aur, it is also un- 
likely that had any extended emission is present. This would 
have been detected. Secondly, electronic transitions of H2 have 
been detected in absorption with the FUSE {Far Ultraviolet 
Spectroscopic Explorer) satellite towards AB Aur , likely arising 
from t he extended envelo pe surrounding the star dRoberge et alJ 
I2OOU iMartin-Zaidi et all [2008 a). The column density of the 
7 = 3 level of H2 measured in the FUSE spectrum is 10^ times 
lower than that derived from the mid-IR observations, implying 
that we clearly do not observe the same regions around the star. 

The source HD 97048 has never been observed with FUSE, 
which precludes any other analysis of the H2 lines than that of 
VISIR observations. We thus analyzed the spectra of HD 97048 
and AB Aur obtained with VLT/UVES. For the two stars we 
detected absorption lines of the CH and CH^ molecules (see 
Fig.|6l), as well as absorption lines of Ca I, Ca II and K I. Due to 
the high inclination angles of the disks with respect to the lines 
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Table 6. Column densities, radial velocities and b parameters of 
CH and CH+ for HD 97048 and AB Aur. 



Star 


log N(CH) 


log N(CH+) 


Vrad 


b 




(cm-2) 


(cm-2) 


(kms-') 


(km s-') 


HD 97048 


13.1+0.4 


12.2±0.3 


14.5±0.2 


1.5+0.4 


AB Aur 


12.3±1.4 


13.1±0.1 


15.0±0.1 


1.8±0.3 



The radial velocities are in the heliocentric rest frame. 



of sight, the observed gas cannot be in the disks, but likely arises 
from the surrounding envelope. We measured the radial veloci- 
ties of these lines in the UVES spectra of the two stars (Table|6]l. 
All these lines in the UVES spectra have the same radial velocity, 
implying that we probed a single gaseous component along each 
line of sight. This is confirmed by the measured b-valnes of CH 
and CH^ lines that are in the range expected for pure thermal 
broadening, as generally seen in diffuse ISM c louds where typ - 
ical Zj-values are about 1.5 up to 3 kms"' (e.g. lGrv et alJl20()2h . 
and which would be much higher if several gaseous components 
were present along the line of sight. In addition, the radial ve- 
locities of the S(l) fines observed for HD 97048 and AB Aur , 
with VISIR and TEXES respectively dMartin-Zaidi et all 120071: 
Bitneret al. 2007) differ by about 10 kms"' from those mea- 
sured for the gas in the UVES spectra. This is an additional clue 
that we probed different regions of the environment of the stars 
in the mid-IR and visible ranges, and that the H2 emission in the 
mid-IR does not come from the envelope, but from the disk. 

Finally, it is important to note here that in their sample 
of Herbig Ae/Be st ars obse rved both with FUSE and UVES, 
iMartin-Zaidi et all (l2009bld) showed that when H2 and CH 
are both observed, their column densities are correlated in 
proportions consistent with those observed in the diffuse in- 
tersteUar me dium (N(CH)/N(H2)~10-\ e.g. [Welty et al. 20(3 
ISheffer et alj |2008). It is important to note that for these stars, 
including AB Aur, the FUSE H2 data implied an interstellar ori- 
gin for the gas or a remn ant of the parent molecul ar cloud and/or 
circumstellar envelope (Martin-Zaidi et al. 2008a). Here, the to- 
tal column densities of H2 derived from the mid-IR observations 
give a N(CH)/N(H2) ratio higher by a factor 10^ to 10^ than the 
interstellar ratio, again suggesting a different location, the H2 de- 
tected in the mid-IR being in the disk, CH being in the envelope. 
Cold molecular hydrogen is very likely present in the envelope 
as probed by FUSE, but does not possess the physical properties 
required to emit at mid-IR wavelengths. 

4.4. Comparison with CO observations 

An important tracer of gas in the inner disks of T Tauri and 
Herbig Ae/Be stars are the fundamental CO ro-vibrational emis- 
sion bands (e.g., Av=l) present at 4.7yum. CO emission at 4.7 
yum traces warm and hot gas at t emperatures rangin g from hun- 
dreds to a few thousands Kelvin (iNaiita et al.ll2003h . For Herbig 
Ae stars, CO has been detected in almost all sources with op- 
tically thick disks (i.e. E(K-L) > 1 Brittain et al. 2007), typi- 
cally exhibiting rotational temperatures b etween 1000-1500 K 
and in a few cases hotter than 2000 K (iBlake & Boogertll2004l: 
Brit tain et al. I l2007h . In Herbig Ae stars, the CO emission prop- 
erties and excitation mechanism appear to depend on the disk 
geo metry deduce d from the SED shape (Meeus et al. 200^ 
Dull emond & Dom inik 2004). The dominant excitation mech- 
anism appears to be LTE for self-shadowed disks, and fluores- 
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Fig. 6. Fits of the CH line at 4300.3 A in the UVES spectra of HD 97048 (left) and AB Aur (right) with the Owens profile fitting 
procedure. Stellar continuum: dashed line; intrinsic Une profile: dotted line; resulting profile convolved with the line-spread function: 
thick line. 



cence for flaring disks (van der Plas et al. 2010, in prep). High- 
excitation CO emission lines (up to v=5-4) have been observed 
only in flared disks (van der Plas et al. 2010, in prep). 

In the context of H2 mid-IR emission detections, we find that 
in all the Herbig Ae/Be stars where H2 mid-IR emission has been 
detected their disks are flared, and CO emission at 4.7/vm has 
been detected too. Note however that the opposite is not true, as 
this and other studies show, not all th e sources with CO 4.1 
emission display H2 mid-IR emission. iBitner et al.l (l2008h show 
that when H2 mid-IR is detected, the gas temperatures are sig- 
nificantly higher than the equilibrium temperatures expected for 
the emitting regions, thus suggesting that the gas temperature is 
higher than the dust and that processes such as UV-X-rays or 
accretion heating may b e important (a conclu sion in agreement 
with the conclusions of ICarmona et al.l I2OO8I when explaining 
non-detections of H2 emission). 

If we e xtend the idea of a n additional so urce of heating to 
collisions jBitner et al.ll2008t ICarmona et al.l l2008) to analyze 
the CO and H2 mid-IR detections together (note that mid-IR 
H2 traces warm gas and CO 4.7 yum hotter gas), we may ex- 
pect that the sources with detected H2 mid-IR emission would 
have CO emission at temperatures higher than the equilibrium 
temperature for the respective emitting regions. In other words, 
we will expect for sources with H2 mid-IR emission to ob- 
serve CO lines at high excitation transitions. This argumen- 
tation works at least for H D 97048, where CO em ission up 
to v-5-4 has bee n detected dvan der Plas et al.l 1200 9*). Indeed, 
Ivan der Plas et al.l (2009) modeled the '^CO emission at 4.7yum 
around HD 97048, observed with the VLT/CRIRES, with a 
homogeneous Keplerian disk. Their modeling of CO emission 
yields an inner radius of 11 AU. Those authors concluded a de- 
pletion of CO in the inner region of the disk, probably due to 
photodissociation. Our previous estimate of the location of the 
H2 mid-IR emission (>5 AU to 35 AU) is compatible with the 
CO inner radius, but CO also traces gas in the disk at much larger 
radii (up to 100 AU). However, this line of reasoning cannot ex- 
plain the case of HD 100546, a star where van der Plas et al. 
(2010, in prep) detected CO emission up to v=5-4, but we do 
not see H2 mid-IR emission with VISIR. But then, there is a 
large gap (r~14 AU) in the dust disk of HD 100546, with very 



little material inside dBouwman et alj2003l;lBenistv et al.l2010l) . 
This is likely different from HD 97048 and AB Aur. 

In addition, iBrittain et alj (l2003h conducted observations 
toward AB Aur at 4.7/im of '^CO (v=l-0) emission using 
the NASA Infrared Telescope Facility (IRTF) and the Keck 
Observatory. The excitation diagram of CO (assuming UTE gas) 
gives a temperature of 70 K for the cool gas coiTesponding to 
the low- J levels, and an excitation temperature of 1540 K for the 
hot gas (high-J levels). This diagram is interpreted by those au- 
thors as hot CO emission coming from the inner rim of the disk 
(~1 AU from the star), while the cool emission originates in the 
outer flared part of the disk (> 8 AU from the star). The excita- 
tion temperature and radial velocity foun d for the cold CO are 
fully consistent with those found by Martin-Zaidi et al.l (l2008ah 
for the cold H2 observed in absorption by FUSE, a s well as for 
the CH observed with UVES. On the other hand, iBitner et"an 
(2007) showed that the H2 emission at 17jum in the disk of 
AB Aur arised around 18 AU from the central star, with an exci- 
tation temperature of about 670 K. Their derived gas temperature 
and distance from the central star fall bet ween the hot and col d 
components seen in the CO observations ( iBrittain et alj|2003l) . 
suggesting that the H2 and CO IR-emission does not come from 
the same region of the disk. 

It is hard to drive definitive conclusions based on very few 
detections. Further searches of high-excitation CO emission us- 
ing high-spectral resolution in sources with H2 mid-IR detec- 
tions (most notably AB Aur) and further searches of H2 mid-IR 
emission in sources with high-excitation CO will be required to 
test the idea whether H2 mid-IR emission and high-excitation 
CO emission are coiTelated. For instance, the main conclusion 
that we can aiTive at is that Herbig Ae stars appear to be in gen- 
eral a uniform group concerning CO fim emission and H2 emis- 
sion in the sense that CO is detected in optically thick disks and 
that H2 mid-IR emission is very weak (i.e. not detected), except 
for those sources with flaring disks where an additional mech- 
anism heats the H2 to a detectable level. The use of full disks 
models (coupling gas and dust; Woitke et al. 2010) could help 
us to understand the possible correlation between CO and H2, 
but this is beyond the scope of this paper. 
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5. Conclusion 

We reported here on a search for the H2 S(l) emission Hne 
at 17.035 jum in the circumstellar environments of five well 
known HAeBes with the high resolution spectroscopic mode of 
VISIR. No source shows evidence for H2 emission at 17.035 
yum. From the 3cr upper limits on the integrated line fluxes, 
we derived limits on column densities and masses of warm 
gas as a function of the temperature. The present work brings 
to 18 the number of non-detections of the H2 S(l) line in 
a global sample of 20 Herbig Ae/Be stars observed with 
VISIR and TEXES dMa rtin-Zaidi et al. ll2007l:lBitner et alj|2007l: 
Carmona et aH 120081: Martin-Zaidi etal.l l2008bt iBitner etall 



2008ir 



Th e det ections of H2 emiss ion at 17.035 fim bv lBitner et alj 
(l2007h and iMartin-Zaidi et al.1 ((2007) show that at least a few 
circumstellar disks have sufficiently high H2 mid-infrared emis- 
sion to be observed from the ground. The most likely explana- 
tion for this is that the optically thin surface layers of the disk 
has Tgas > Tdust and th at the gas-to-d ust ratio i s higher than the 
cano nical ratio of 100 dCarmona et al . 2008; Martin-Zaidi et alj 
l2007h . Indeed, in the surface layers of the disk, low densities or 
dust settling and coagulation may conduct to a spatial decou- 
pling between the gas and the dust. Photoelectric heating can 
thus play a significant role in the gas heating process and the 
physical con ditions may rapidly differ from the LTE ones, i.e ., 
Tgas > Tdust dKamp & Dullemondlllool iJonkheid et a"l]|2007h . 
On the other hand, UV and X-ray heating can be responsible for 
the excitation of the observed gas and c an heat the gas to temper- 
atures significantly ho t ter than the dust ( Nomura & Millaij[2005l: 
lGlassgoldetani2007t lErcolano etaDl2008h . In any case, one 
would need to observe the lower H2 7-levels (i.e. 7 = and 
J - 1) to definitively constrain the kinetic temperature of the 
gas, and better understand the excitation mechanisms responsi- 
ble of the mid-lR emission. 

As a second step, we also performed a statistical analysis 
of the whole sample of Herbig Ae/Be stars observed at 17.035 
fj.m with VISIR as well as with TEXES. This analysis allowed 
us to identify a population of stars, including HD 97048 and 
marginally AB Aur, with properties that depart from the bulk 
of our sample. This raises numerous questions about the origin 
of the detected gas and the status of HD 97048 and AB Aur. 
From our UVES observations we clearly demonstrated that the 
observed mid-lR H2 emission does not come from the envelope, 
but from the disk. 

Are HD 97048 and AB Aur peculiar stars (and why)? Due to 
the similarities of these two stars (Te/f, age, mass, disk size...), 
one would expect that the physical conditions of their circum- 
stellar gas are typical of a particular (short) stage of evolution 
of the disks. However, to confirm this assumption, we would 
need to observe other similar HAeBes. Unfortunately, no other 
nearby Herbig Ae/Be star observable with the existing instru- 
ments presents the same observational properties. A global diag- 
nostic of the gaseous content of the disks is thus now required. 
To better constrain the physical and chemical properties of the 
gas, multi-wavelengths observations and a deep modeling would 
be very useful. In this context, the Herschel satellite will be very 
helpful because it will allow us to constrain the thermodynamics 
in young disks by observing the gas and the dust simultaneously. 
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